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INTRODUCTION 


During  this  quarter  cur  efforts  have  been  directed  toward 
integrating  the  propellant  response  model  with  the  experimental  results 
in  a  more  basic  way. 

A  gas  dynamic  study  has  shown  tba t  the  finite  amplitude  traveling 
wave  can  be  sustained  in  a  combustion  chamber  hy  a  modest  input  of 
energy — 1 %  to  3%  of  the  combustion  energy;  the  source  of  this  energy 
might  be  easily  explained  by  partially  completed  reactions  at  the  sur¬ 
face.  However,  the  fact  that  there  are  certain  limits  to  the  frequency 
at  which  instability  occur*  suggests  that  gas  phase  perturbations  respond 
tc  reactions  with  a  relatively  long  relaxation  time.  The 
analytical  study  has  revealed  that  a  large  response  is  associated  with 
a  significant  fractional  heat  release  <for  composite  propellants) 
occurring  in  the  condensed  phase  or  surface-coupled  reactions.  In  the 
experimental  study,  attempts  are  being  made  to  identify  the  nature  of 
solid  phase  reactions  and  to  differentiate  between  condensed  phase  and 
gas  phase  reactions.  Adiabatic  self-heating  studies  and  diiierentifi 
thermal  analysis  studies  are  being  performed  at  a  range  oi  pressures. 

Preliminary  experimental  results  suggest  that  between  10%  and  20% 
of  the  heat  release  oceurs  in  the  solid  phase  reactions  of  certain  A P~ 
based  propellants;  the  condensed  phase  heat  release  for  KP-based 
propellants  appears  to  be  much  less. 

U3ing  activation  energies  obtained  in  experimental  studies  prev¬ 
iously  reported  in  this  program,  it  has  proved  possible  to  theoretically 
predict  the  frequency  band  for  propellant  response.  The  magnitude  of 
this  response  is  dependent  on  the  heat  release  in  the  solid  (condensed) 
phase.  The  frequency  band  predicted  agrees  closely  with  the  preferred 
frequencies  experimentally  observed.  This  appears  to  represent  a  sig¬ 
nificant  advance  In  theoretical  prediction  of  unstable  combustion 
characteristics  from  ballistic  data  (burn.'ng  rate)  and  physical  chemical 
data  (activation  energy). 
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THEORETICAL  INTERPRETATION  OF 
EXPERIMENTAL  RESULTS 


Traveling  Lave  Phenomena 

During  the  past  quarter  attention  has  been  focused  on  the  behavior 
of  traveling  waves  in  the  combustion  chamber,  and  especially  upon 
further  interpretation  of  the  available  data.  Typical  experimental 
behavior  is  shown  in  Figure  1  for  PBAN  319  propellant  in  a  slab  motor 
operating  at  a  mean  pressure  of  1260  psia.  The  important  points  to  be 
noted  are  that  the  head-end  and  aft-end  pressure  pulses  are  nearly 
equal  in  magnitude  and  approximately  twice  as  large  as  the  pulses 
measured  at  the  one-quarter  and  three-quarter  length  stations.  The 
aft-end  pulse  is  slightly  smaller  than  the  head-end  pulse,  presumably 
because  of  some  energy  loss  upon  reflection  from  the  open  nozzle. 

From  a  phenomenological  point  of  view,  the  behavior  illustrated 
by  the  pressure  traces  of  Figure  1  can  be  explained  by  the  presence 
of  a  constant  strength  shock-expansion  process  which  is  traveling  back 
and  forth  in  the  motor.  The  passage  of  the  shock  wave  past  any  point 
induces  a  particle  velocity  behind  it  in  the  direction  of  travel  of 
the  shock,  as  well  as  pressure  and  temperature  jumps.  In  order  to 
satisfy  the  continuity  equation,  an  expansion  process  which  reduces 
the  induced  velocity  to  zero  must  form  behind  the  shock  wave,  as  shown 
schematically  in  Figure  2.  Since  the  expansion  process  is  isentropic, 
the  local  velocity  of  the  expansion  field  will  be  the  local  speed  of 
sound  which  decreases  with  increasing  distance  behind  the  shock  wave. 
Thus  the  extent  or  me  expansion  process  will  lengthen  as  wave  travel 
proceeds  down  the  chamber. 

When  the  shock  wave  reflects  from  the  end  of  the  chamber,  the 
measured  perturbation  amplitude  doubles  because  the  shock  wave  will 
maintain  the  same  pressure  change  across  itself,  whereas  the  pressure 
in  front  of  the  reflected  shock  is  the  pressure  which  was  behind  the 
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FIG.  1  TRAVELIN 


WAVE  INSTABILITY  IN  A  SLAB  BURNER 
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-..cident  shock.  A  correspondingly  strong  expansion  process  will  follow 
the  reflected  shock  back  down  the  passage  toward  the  oncoming  expansion 
process  which  was  following  the  incident  shock.  This  incoming  expansion 
process  reduces  the  strength  of  both  the  reflected  shock  and  its  trailing 
expansion  process  down  to  the  strengths  of  the  original  process  and  the 
chain  of  events  is  repeated  in  the  opposite  direction.  This  complicated 
turn-around  process  occurs  in  a  very  short  physical  distance ,  making 
the  details  very  difficult  to  obtain  experimentally. 

This  physical  reasoning  is  borne  out  by  further  examination  of 
Figure  1.  Not  only  does  pressure-doubling  occur  at  the  ends,  but  the 
wave  velocity  along  the  chamber  is  nearly  constant,  since  a  comparison 
of  the  data  at  the  quarter-length  point  with  that  at  the  head-end  shows 
that  about  one-quarter  of  the  total  wave  travel  time  is  used  over  the 
first  quarter  of  the  motor.  A  constant  wave  velocity  implies  a  constant 
wave  strength ;  within  the  accuracy  of  the  traces ,  the  strengths  at  the 
one-quarter  and  three-quarter  points  are  the  same. 

The  final  strength  (and  therefore  velocity)  of  the  wave  will  be 
determined  by  the  balance  between  the  energy  input  process  and  the 
dissipation  process  at  the  walls  and  by  the  coupling  between  the  passage 
of  the  wave  and  the  primary  combustion  process.  This  latter  coupling 
is  discussed  in  detail  in  a  later  section.  A  possible  mechanism  for  the 
local  input  of  energy  into  the  traveling  wave  is  the  heat  release  from 
gas-phase  reactions  which  are  induced  by  the  temperature  jump  across 
the  shock  wave. 

If  the  traveling  wave  is  considered  as  a  Chapman-Jouget  detonation, 
the  pressure  perturbation  is  related  to  the  Mach  number  and  the  pressure 
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and  the  heat  release  required  to  drive  the  wave  is 
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For  the  data  of  Figure  1,  Ap/pj  *=  0.093  giving  K  *=  1.08  and  Q/c  Tx  = 

1.25  x  10  .  Since  c  T,  is  of  the  ordc’'  of  the  available  energy  per 

P 

unit  mass  of  the  propellant,  it  is  cleer  that  only  a  small  fraction  of 
this  available  energy  is  required  to  drive  the  shock  wave.  It  must 
be  noted  that  the  above  calculation  does  not  account  for  energy  lost 
in  dissipation  at  the  walls,  so  that  the  actual  Q  extracted  from  the 
stream  is  larger  than  the  value  calculated  above.  Since  the  dissipation 
might  be  expected  to  be  less  in  a  tubular  burner  than  in  a  slab  burner, 
energy  should  be  available  to  drive  a  stronger  wave  in  a  tubular  burner. 
That  this  is  indeed  the  case  is  shown  belcw. 

To  summarize,  traveling  wave  data  obtained  in  a  slab  burner  supports 
the  concept  of  the  instability  being  a  constant  strength  detonation  wave 
supported  by  an  energy  input  from  gas  phase  reactions  which  is  small 
compared  to  the  total  energy  available  in  the  propellant.  More  important 
is  the  coupling  between  the  wave  process  and  the  main  combustion  process 
which  is  discussed  in  detail  below. 
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Tubular  Motor  Results 

Head-end  pressure  measurements  obtained  in  tubular  motors  with 
five  different  propellants  are  shown  in  Figure  3.  Two  traces  are  shown 
for  each  propellant,  the  lower  one  of  which  har  been  filtered  to  remove 
the  organ  pipe  oscillation  in  the  small  cylindrical  chamber  ahead  of 
the  transducer.  Results  are  given  in  Table  I. 

The  wave  amplitude  (see  Table  I),  computed  as  half  the  head-end 
pressure  pulse  for  reasons  discussed  above,  is  relatively  constant  for 
all  the  ■propellants  considered,  being  about  Ap/p  =  0.17.  The  variations 
which  do  occur  undoubtedly  depend  more  upon  the  ratio  of  the  mean  pressure 
to  the  threshold  pressure  for  instability  than  upon  compositional  factors 
in  these  similar  AP-based  propellants.  Using  this  value  for  Ap/p, 
and  assuming  that  p  =  p  -  (Ap/2),  equation  (1)  gives  M  ~  1.2  for  the 
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FIG.  3  HEAD-END  PRESSURE  TRANSIENTS  DURING  TRAVELING  WAVE 
INSTABILITY  IN  A  TUBULAR  BURNER 
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FIG.  3  (Continued) 
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TABLE  I 


TUBULAR  MOTOR  INSTABILITY  DATA 


Propellant 

Composition 

Mot  or 

Length 

(inches) 

Wave 

Amplitude 

(Ap/p) 

Wave 

Frequency 

(cps) 

PBAN  103 

24%  ground* 

,  I  NIL  CIO. 

56%  unground  1  4  4 

20%  PBAN 

40 

0.185 

600 

PBAN  104 

22.5%  ground*^ 

f  nh.  CIO. 

56%  unground  j  4  4 

20%  PBAN 

1.5%  LiF 

82 

0.169 

244 

PBAN  244 

27%  ground*  'I 

52%  unground  J  N^tcl’04 

20%  PBAN 

1%  LiF 

40 

0.163 

500 

PBAN  284 

20%  -  20p 

20%  -  600p.  InF^CIO* 

39.5%  ungrounrij 

20%  PBAN 

0.5%  SrC03 

40 

0.138 

513 

PBAN  319 

I 

1 _ _ 

22%  -20^,  > 

22%  -  600y,  \  N^C104 

34%  un ground  J 

17%  PBAN 

£%  A1 

82 

0.184 

250 

♦Average  particle  diameter  of  ity. 
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wave  Kach  number  and  equation  (2)  gives  Q/c  /T,  =  6.3  x  10  .  Even 

P 

accounting  for  dissipation,  it  conservatively  appears  that  less  than 
2  percent  of  the  propellant  energy  is  needed  for  direct  input  to  the 
wave.  Notice  that  the  dissipation  is  much  smaller  in  the  tubular  burner 
than  in  the  slab  burner. 


INTERACTION  OF  THE  TRAVELING  WAVE  AND  THE 
COMBUSTION  MECHANISM 


As  has  been  well-documented  in  this  program  and  others,  axial -mode 
combustion  instability  arises  when  conditions  are  such  that  a  pressure 
disturbance  does  not  decay,  but  instead  develops  into  a  traveling 
wave  reflecting  back  and  forth  between  the  ends  of  the  chamber.  The 
preceding  discussion  shows  that  this  wave  is  actually  a  weak  shock  wave 
its  strength  depends  to  some  degree  on  the  propellant  formulation  and 
chamber  conditions,  but  is  almost  independent  of  axial  position  in 
the  chamber.  Usually,  the  Mach  number  of  this  traveling  shock  wave 
is  approximately  1.2. 


The  energy  feedback  from  the  reaction  zone  required  to  maintain 
a  shock  of  this  strength  represents  only  about  2  percent  of  the  total 
rate  of  energy  release  in  the  typical  propellant  combustion  process, 
This  observation  leads  to  a  very  important  conclusion;  it  is  unlikely 
that  the  appearance  of  axial-mode  instability  is  governed  primarily 
by  the  amount  energy  feedback  from  the  combustion  process  to  the 
traveling  wave,  since  so  little  energy  is  needed  to  support  the  wave. 
Of  course,  tfc «  energy  flux  magnitude  does  play  a  significant  role, 
particularly  in  a  motor  with  a  high  damping  factor,  or  high  losses. 
However,  it  is  likely  that  a  more  important  factor  in  determining  the 
axial-mode  stability  characteristics  of  a  solid  rocket  motor  is  the 
phase  relationship  between  the  energy  feedback  at  any  given  point 
on  the  propellant,  and  the  passage  of  the  wave  over  that  point.  In 
other  words,  if  the  energy  feedback  (or  mass-flux  perturbation  at  the 
surface,  which  is  essentially  equivalent)  is  locally  in  phase  with  the 
passing  shock  wave  over  most  of  the  burning  surface,  the  shock  wave 
will  be  reinforced  and  instability  will  result.  If  the  two  events 
are  locally  out  of  phase  over  most  of  tbs  longth  of  the  grain,  the 
shock  wave  will  decay  and  stable  combustion  will  be  achieved.  This 
phase  relationship  is  determined  by  the  combustion  chamber  length  and 
by  the  frequency  response  characteristics  of  the  combustion  mechanism. 


Local  Frequency  of  the  Pressure  Pulse  Caused  by  the  Passing  Shock  Wave 


As  the  shock  wave  travels  axially  back  and  forth  in  the  combustion 
chamber,  it  creates  an  oscillatory  pressure  disturbance  at  every  point 
on  the  propellant  surface.  The  frequency  of  this  disturbance  at  any 
given  axial  position  corresponds  to  the  number  of  times  per  second 
chat  the  wave  passes  this  position.  In  general,  there  are  two  distinct 
frequencies  associated  with  any  given  axial  position.  One  of  these 
corresponds  to  the  time  required  for  a  shock  wave  which  is  traveling 
toward  the  right,  to  be  reflected  from  the  right  end  of  the  chamber  and 
return  to  the  reference  point.  The  other  frequency  corresponds  to  the 
time  for  a  wave  traveling  toward  the  left  end  to  be  reflected  and  return 
to  the  some  point.  At  the  ends  of  the  chamber  only  one  of  these  frequencies 
pertains  (i.e.,  the  other  is  infinite);  at  the  center  they  are  equal, 
of  course.  From  a  brief  consideration  of  this  phenomenon  one  can  readily 
show  that  the  pulse  frequency  at  the  center  of  the  chamber  1b  exactly 
twice  that  at  the  ends.  At  intermediate  stations  the  propellant  iB 
subjected  to  a  train  of  pulse  pairs.  For  example,  at  one-quarter  and 
three-quarter  positions  the  separation  between  pairs  corresponds  to  a 
recurrence  frequency  of  two-thirds  of  that  at  the  center,  while  the 
spacing  between  two  pulses  in  a  pair  corresponds  to  a  frequency  of 
twice  that  at  the  center.  The  band  of  frequencies  encountered  is 
dependent  on  the  motor  length  and  the  motor  surface  is  dependent  on 
the  net  driving  and  damping  of  the  wave  train. 


Pulse  frequencies  encountered  in  experiments  performed  during 
this  program1  are  shown  in  Table  II. 


For  the  ?3-,  40- ,  and  82-inch  motors,  frequencies  were  measured 
at  the  head-end  pressure  transducer.  As  was  shown  in  earlier  discussion 
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wave  moving  at  M  =*  1.2  back  and  forth  in  the  chamber.  Ordinarily, 
pulses  introduced  in  the  15-inch  motor  were  found  tc  decay,  i.e.,  with 
most  propellantB  this  motor  could  not  be  driven  tc  unstable  combustion 
The  frequencies  shown  in  Table  II  for  this  motor  are  the  calculated 


1Capener,  E.  L  .  R.  J.  Kier,  L.  A.  Dickinson,  and  G.  A.  Marxman, 
Response  of  a  Burning  Propellant  Surface  to  Erosive  Transients,  Final 
Scientific  Report,  A.F.  Office  of  Scientific  Research  Contract 
AF  49(638)-1507 ,  March  15.  1966. 


TABLE  II 


FREQUENCY  OF  PRESSURE  PULSE,  CPS 


Length  of 
Chamber,  Inches 


2480 


characteristic  frequencies  based  on  a  traveling  shock  wave  with  VL  -  l.‘. 
Also,  it  is  important  to  note  that  in  sco>^  cases  the  t>2”  motor  tended 
to  go  unstable  in  a  double  mode,  i.e.  with  traveling  waves  in  two 
4l"  sections  reflecting  against  each  other  at  the  center.  In  this 
case  the  frequencies  for  the  82”  motor  are  essentially  the  same  as  for 
the  40"  motor,  where  the  center  of  the  longer  motor  acts  as  the  “'end” 
of  two  shorter  ones. 


Theoretical  Trequency  Response  of  the  Propellant  Combustion  Mechanism 


As  a  f.'rst  step  toward  explaining  our  experimentally  determined 
limits  of  axial-mode  combustion  instability,1  it  is  essential  to  com¬ 
pare  the  observed  frequency  characteristics,  as  reported  in  Table  I, 
with  the  frequency  response  characteristics  of  the  propellant  combustion 
mechanism.  The  latter  has  been  predicted  by  the  theoretical  combustion 
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3Capener,  E.  L. ,  L.  A.  Dickinson,  and  G.  A.  Karxman,  Response  of  a 
Burning  Propellant  Surface  to  Erosive  Transients,  First  Quarterly 
Report,  Air  Force  Office  of  Scientific  Research,  Contract  No. 

AF  49(638) -1665,  April  1966. 

3Capener,  E.  L. ,  I,.  A.  Dickinson,  and  G.  A.  Marxman,  Response  of  a 
Burning  Propellant  Surface  to  Ero3ivo  Transients,  Second  Quarterly 
Report,  Air  Force  Office  of  Scientific  Research,  Contract  No. 

AF  49(638'  1665,  July  1966. 
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1/4  L 

3/4  L 

1660 

4960 

1070 

3212 

620 

1860 
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observed  instability  characteristics  in  terms  of  the  analytical  predic¬ 
tions  affords  a  significant  test  of  the  underlying  theoretical  concepts. 

As  a  specific  example,  Figure  4  shows  the  theoretical  frequency- 
response  characteristic  of  a  composite  propeller*  in  which  about  5 
percent  of  the  total  heat  release  occurs  in  surface-coupled  reactions, 
and  the  rest  occurs  in  the  usual  gas-phase  reactions.  This  figure  is 
reproduced  from  an  earlier  report,3  wherein  its  basis  is  described.  The 
coordinates  are  composed  of  groupings  of  thermo-chemical  parameters  of 
the  propellant.  Thus,  any  propellant,  having  certain  values  of  the 
activation  energies,  surface  temperature,  etc.,  can  be  represented 
as  a  point  of  Fig.  1.  The  solid  curves  show  the  maximum  ratio  of 
burning  rate  oscillation  amplitude  to  pressure  disturbance  amplitude 
that  can  occur  at  any  given  point.  Only  a  certain  frequency  of  pressure 
disturbance  (shown  by  the  dashed  curves)  can  induce  that  response.  Any 
other  frequency  of  pressure  oscillation  causes  a  smaller-amplitude 
burning-rate  oscillation. 

The  first  thing  to  note  about  Fig.  4  is  that  the  frequency  at 
which  the  propellant  responds  most  strongly  depends  almost  entirely  on 
the  value  of  A  (i.e.  the  dashed  curves  are  nearly  horizontal) .  This 
is  even  more  true  at  higher  values  of  9  (greater  percentages  of  surface- 
coupled  heat  release)  than  that  associated  with  Fig.  1.  Thus,  the 
resonant  frequency  of  the  combustion  mechanism,  f  =  u/2tt,  depends  on 
the  burning  rate  r,  the  thermal  diffusivlty  of  the  solid  K,  and  also 
on  the  activation  energy  of  surface  decomposition  E  as  well  as  the 
surface  temperature . 

The  other  important  point  shown  by  Fig.  4  is  that  the  maximum 
amplitude  of  the  burning  rate  oscillation,  as  represented  by  the  amplifi¬ 
cation  factor  't/p.  depends  almost  entirely  on  the  value  of  ™  atmru^i  atari 
with  that  propellant.  (This  is  not  true  at  low  values  of  A,  but  few 
actual  propellants  fall  in  that  range.)  Thus  the  maximum  possible 
amplification  factor  for  a  given  propellant  is  determined  largely  by 
the  activation  energy  of  the  gas-phase  flame  Ef ,  the  gas-phase  flame 
temperature  T^,  and  the  effective  overall  order  n  of  the  gas-phase 
react  ions . 
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However,  there  is  a  third  important  point  that  cannot  be  deduced 

fraai  Fig.  4  alone;  the  amplitude  of  the  burning  rate  response  is  also 

strongly  influenced  by  ©^.  It  can  be  shown3  that  as  ©s  increases 

(increasing  percentage  of  surface-coupled  heat  release)  the  solid 

curves  move  toward  the  left,  and  the  "bump”  in  these  curves  in  the 

region  1  <  A  <  3  is  eliminated.  (This  occurs  rapidly;  the  bump  is 

gone  at  ©^  «  1.0).  This  means  that  of  two  propellants  with  similar 

values  of  the  theraocbemical  parameters  (and,  therefore,  the  values 

of  A  and  a) ,  the  one  with  more  surface -coupled  heat  release  (greater 

0  )  will  have  a  higher  amplitude  burning-rate  oscillation  for  a  given 
© 

pressure  oscillation.  Inasmuch  as  most  propellants  probably  have 
similar  values  ol’  a,  this  effect  may  be  a  key  factor  it-  separating 
those  that  burn  stably  from  those  that  do  not. 


To  put  these  theoretical  conclusions  in  a  more  quantitative  form 

it  is  necessary  to  assign  approximate  values  to  the  thermochemical 

parameters  that  characterize  the  propellants.  For  most  composite 

propellants  it  is  reasonable  to  assume  that  the  activation  energy 

for  surface  decomposition  is  E  w  30  kcal.  The  surface  temperature 

T  is  approximately  800°K,  and  the  gas-phase  flame  temperature  T.  2500°K 
w  I 

The  activation  energy  for  the  gas-phase  flame  usually  is  E^  m  20  kcal 
and  the  overall  effective  order  of  that  reaction  n  is  about  2.  The 
parameters  for  a  very  wide  range  of  propellants,  including  those  that 
have  been  investigated  experimentally  in  this  program,  almost  certainly 
fall  within  20%  of  these  figures.  Thus,  a  typical  lomposite  propellant 
is  characterized  in  Fig.  4  by  the  coordinates  A  ~  10,  1/a  ~  1.  It 
follows  from  the  theory  (Fig.  4)  that  the  resonant  frequency  f  of  the 
combustion  mechanism  for  such  propellants  corresponds  to  Ku/r3  ~  9, 
or  f  ~  9t9/2ttK.  It  also  can  be  shown  that  tbs  maximum  deviation  from 
this  frequency  for  which  there  is  a  significant  burning-rate  amplifica¬ 
tion  is  about: 


(This  frequency  range  follows  directly  from  an  examination  of  the  mathe¬ 
matical  analysis  and  may  be  deduced  by  referring  to  the  mathematically 
similar  treatment  by  Denison  and  Baum.4)  A  typical  thermal  diffusivity 
for  solid  propellants  is  K  m  2.3  x  lO-"*  inVsec.  Accordingly,  the 
approximate  frequency  range  in  which  the  combustion  mechanism  tends  to 
amplify  a  pressure  oscillation  is  as  indicated  in  Table  III.  The 
frequency  given  in  the  table  is  the  approximate  resonance  frequency,  or 
point  of  maximum  amplitude,  and  the  tolerance  figures  indicate  the 
approximate  range  in  which  a  degree  of  amplification  is  present. 


TABLE  III 

APPROXIMATE  RESONANT-FREQUENCY  RANGE  FOR  A 
TYPICAL  COMPOSITE  PROPELLANT  ACCORDING  TO  THEORETICAL  COMBUSTION  MODEL 


r  (in/sec) 


£  (cps) 

~  60  ±  15 
~  250  ±  80 
~  560  ±  130 
~  1000  ±  210 
~  1550  ±  340 


It  is  reasonable  to  assume  that  axial-mode  combustion  Instability 
may  arise  in  a  solid  rocket  chamber  whenever  an  appreciable  portion 
of  the  grain  length  is  subjected  to  pressure  oscillations  near  the 
resonance  frequency  cf  the  combustion  mechanism.  Cn  this  basis  a 
comparison  of  Tables  II  and  III  afford  some  interesting  conclusions. 

In  a  15-inch  chamber,  only  a  very  high  burning  rate  propellant 
should  be  susceptible  to  axial-mode  instability.  With  a  0.5-inch/sec 
burning  rate,  only  +he  portion  of  the  propellant  near  the  ends  of  the 
chamber  can  experience  pressure  oscillations  in  the  resonance  range. 
With  lower  burning  rates,  no  part  of  the  comb'  stion  surface  is  in 
resonance  with  the  waves.  Thus,  with  propellants  in  the  normal  burning 
race  range  the  theory  predicts  that  combustion  instability  should  be 
difficult  to  initiate  (or  rarely  experienced)  '5-incb  chamber. 

This  agrees  with  our  experimental  observations. 

4 M  R.  Denison  and  E.  Baum,  "A  Simplified  Model  of  Unstable  Burning  in 
Solid  Propellants,”  ARS  Journal,  31,  1112  (1961). 


— — 


Instability  is  much  more  likely  in  the  23-inch  and  40-inch 

ehambe  s.  In  the  former,  propellants  with  burning  rates  of  0.35  to 

0.45  inch/aec  should  be  particularly  susceptible  to  instability.  The 

corresponding  range  in  the  40-inch  chamber  is  about  C.2  to  0.4 

inch/sec.  The  stability  characteristics  of  a  number  of  propellants 

have  been  thoroughly  documented  experimentally  in  a  40-inch  chamber 

during  this  program.1  The  results  are  shown  in  Fig.  5.  It  is  apparent 

that  the  burning  rates  of  the  propellants  that  encountered  instability 

fell  in  the  range  0.1  to  0.3  inch/sec.  Though  this  range  is  somewhat 

lower  than  that  predicted  theoretically,  the  agreement  is  well  within 

the  tolerances  imposed  by  uncertainties  in  evaluating  the  thermochemical 

constants  of  the  propellant.  The  theoretically  predicted  resonance- 

frequency  bounds  offer  a  plausible  explanation  for  the  fact  that  the 

high  burning  rate  AP/KP  \ .opellaut  was  stable ,  whereas  the  other  one 

was  not.  They  may  also  explain  why  the  KP ,  LiP,  and  AN  propellants  were 

stable,  as  all  these  have  resonance  frequencies  well  outside  the  wave 

frequencies  typical  of  a  40-inch  motor.  However,  there  is  also  another 

factor  that  may  be  significant  here-- the  value  of  6  which  determines 

8 

the  amplitude  of  the  burning  rate  response.  This  factor  certainly 
appears  to  be  important  relative  to  the  opposite  behavior  of  the  AN/KP 
and  AN/AP  propellants,  which  chould  have  similar  resonance  frequencies. 
We  shall  consider  this  point  further  in  later  discussion. 

In  the  82-inch  motor  it  was  found  experimentally  that  the  shock 
wave  traveled  from  end  to  end  of  the  motor,  with  the  frequencies  shown 
in  Table  II,  only  during  the  first  100  msec  after  the  initiating 
pulse.  In  every  case  the  motor  failed  to  sustain  such  a  mode  beyond 
this  initial  period.  Instead,  &  transition  to  a  double  mode,  or 
essentially  the  first  harmonic,  always  occurred.  In  this  case  there 
were  two  traveling  shock  waves  in  the  motor,  each  traveling  through 
a  41-inch  section  and  reflecting  off  the  other  wave  at  the  center. 

Thus,  the  82-ineh  motor  always  undergoes  a  transition  to  an  axial 
instability  mode  that  exhibits  the  same  frequency  range  as  that  shown 
for  a  40-inch  motor  in  Table  II.  This  is  precisely  what  one  would 
expect  for  this  propellant  (PBAN  103,  r  ~  0.25-0.3  inch/sec)  from 
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LINEAR  BURNING  RATE  -  in./«c 


>»AN/AP 


AP  STABILITY  BOUND 


I »-»OU3-3B 


INFLUENCE  OF  BURNING  RATE  AND  COMPOSITION  ON 
FINITE  AMPLITUDE  TRAVELING  WAVE  INSTABILITY; 
SOLID  LINE  STABLE  REGIME,  DOTTED  LINE  UNSTABLE 
REGIME  FOR  5"!NCH  x  40-  INCH  MOTOR 


the  theoretical  analysis  of  Table  III.  The  fundamental  mode  of  the  82- 
inch  motor  corresponds  to  a  frciuency  range  (Table  II)  that  is  somewhat 
below  the  resonanc'  frequency-response  range  of  the  propellant.  However, 
the  double-wave  mode,  whici  corresponds  to  the  frequencies  of  a  40-itich 
motor  (Table  I),  is  right  in  the  resonance  range,  as  has  already  been 
shown.  Thus  the  transition  to  this  mode  in  the  longer  motor  is  entirely 
consistent  w^th  the  theoretical  predictions. 

The  Influence  of  Surface -Coupled  Reactions 

The  foregoing  discussion  deals  only  with  the  resonance-frequency 

aspect  of  the  combustion  model.  To  fully  explain  the  data  it  is 

necessary  to  consider  the  remaining  theoretical  parameter  0g,  which 

is  a  measure  of  the  surface-coupled  heat  elease.  In  general,  when 

0  is  large  for  composite  propellants  (e.g.  10-20%  of  total  energy 
s 

release  in  surface-coupled  reactions)  ,  the  amplitude  of  the  burning 
rate  response  in  the  resonance  frequency  range  is  la>%  e.  When  0g 
is  small  (10%  or  less  heat  release  near  the  surface)  the  amplitude 
is  greatly  reduced,  according  to  the  theory.  This  affords  a  possible 
explanation  for  the  behavior  of  the  AN/KP  and  AN/'AP  propellants  in 
Fig.  5.  It  has  been  shown  during  this  program  that  AP  decomposition 
and  combustion  is  associated  with  a  high  surface  exotherm,  whereas  the 
surface-coupled  energy  release  is  believed  to  be  lower  with  KP.  Thus, 
altoough  both  the  AN/AP  and  the  AS/KP  propellants  are  in  the_x  resonant 
frequency  range  in  the  40-lnch  motor,  the  former  probably  has  **  much 
stronger  response.  Consequently  the  normal  damping  phenomena  in  the 
motor  cannot  prevent  the  AP  composition  from  going  axially  unstable, 
but  they  can  do  so  with  the  KP  propellant.  In  terms  of  Fig.  5,  this 
means  that  the  stability  bound  for  the  AN/KP  propellant  is  approximately 
parallel  to  the  AP  bound  shown,  but  has  a  lower  intercept  on  the 
burning  rate  coordinates.  This  was  derived  from  the  theory  previously  1 


SUMMARY 


The  theoretical  combustion  model  affords  a  plausible,  consistent, 
and  cfmprehensive ,  though  not  proven,  explanation  for  the  experimentally 
observed  axial-raode  instability  phenomena.  Two  factors  enter:  (1)  the 
resonance  frequency  range  of  the  combustion  mechanism,  and  (2)  the 
amplitude  of  the  response  in  the  range.  The  former  depends  primarily 
on  the  burning  rate  and  also  on  thermochemical  characteristics  of  the 
decomposition  process  (activation  energy,  surface  temperature).  It 
can  be  evaluated  quantitatively  with  reasonable  accuracy.  The  amplitude 
response  depends  primarily  on  the  degree  of  surface-coupled  heat 
release,  and  also  on  thermochemical  parameters  of  the  flame.  As  yet 
we  do  not  know  how  to  measure  or  accurately  predict  the  surface-coupled 
reactions  associated  with  a  given  propellant.  Therefore,  the  amplitude 
factor  is  difficult  to  predict  a  priori ,  although  experiments  such  as 
those  reported  in  Fig.  5  are  in  a  sense  a  measure  of  It.  Thi3  is 
certainly  an  area  that  will  require  much  more  attention,  both  experi¬ 
mentally  and  theoretically. 
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